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We derive an effective Maxwell-London equation for entangled polymer complex under the topo¬ 
logical constraint, borrowing the theoretical framework from the topological field theory. We find 
that the transverse current flux of a test polymer chain, surrounded with the entangled chains, 
decays exponentially from its average position with finite penetration depth, which is analogous to 
the magnetic-field decay in a superconductor (SC). Like the mass acquirement of photons in SC is 
the origin of the magnetic-field decay, the polymer earns uncrossable intersections along the chain 
due to the preserved linking number, which restricts the deviation of the transverse polymer current 
in the normal direction. Interestingly, this picture is well incorporated within the most successful 
phenomenological theory of the so called tube model, of which researchers have long pursued its 
microscopic origin. The correspondence of our equation of motion to the tube model claims that 
the confining tube potential is a consequence of the topological constraint (linking number). The 
tube radius is attributed to the decay length, and increasing the retracting force at intersections or 
increasing the number of intersections (linking number), the tube becomes narrow and tighter. It 
further shows that the probability of the tube leakage decays exponentially with the decay length 
of the tube radius. 


Introduction: It is an intriguing nature of the physics that 
more than two completely different systems are described 
within the same mathematical framework. The Chern- 
Simons theory which was first introduced in string theory 
by Witten at I980’s [T], turns out to be applicable to the 
cutting edge problems in condensed matter, such as frac¬ 
tional quantum Hall effect [2] and topological insulator 
and superconductor [5]. Since Edwards first introduced 
the topological constraint [3], i.e. linking number, to the 
partition function of the polymer melts, researchers no¬ 
ticed the potential usefulness to apply the Chern-Simons 
theory to polymer problems [SHE]- In those works, the 
knotted electron’s movement is translated to the inter¬ 
section (entanglement) of polymers. In spite of its math¬ 
ematical exactness, these approaches were not very suc¬ 
cessful in polymer entanglement unlike its huge achieve¬ 
ment in other fields of physics 0110]. Main reasons are 
that the model is still abstract, and the information of mi¬ 
crostates is not feasible under current experimental tech¬ 
niques 0[ini. 

Surprisingly, long polymer melts, which is one of the most 
complicated system in polymer physics, is well described 
by a single body interpretation, i.e. tube model ffn- 
[ 13 ]. In addition to many rheological experimental sup¬ 
ports na EHE], recent experiments and simulations 
verified the validity of the model finding the linear con¬ 
finement potential around semi-fiexible polymers [IIHIS]- 
The strength of this intuitive model lies in its universal¬ 
ity spanning over the intermediate and long time scale 
regardless of the detailed chemistry or the structure of 
the polymer, as long as the polymer is long enough so 
that we can treat it as a coarse grained random chain 
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with the stepsize of random work much exceeding the 
local length scale of the polymer. Since this phenomeno¬ 
logical model earned a great success, researchers have 
sought its microscopic origin j^OHll]. Recent theoretical 
achievements in this effort have focused on the Langevin 
approach of statistical mechanics [331 HI] • 

These two approaches, Chern-Simons theory and the 
tube model, look completely different. However, we be¬ 
lieve that the phenomenological theory should be under¬ 
stood from the microscopic theory. The goal of the cur¬ 
rent work is to pursue this connection. 

As Edwards claimed in his pioneering work, polymer con¬ 
formation can be translated as a classical particle mo¬ 
tion [IlHS]. In this analogy, time corresponds to the 
displacement along the chain, and the trajectory of the 
particle motion is associated with the contour of the poly¬ 
mer. Within this argument, the diffusive motion of a 
classical particle can be interpreted as a Gaussian distri¬ 
bution of a random polymer chain. To include the entan¬ 
glement effect, he plugged the topological constraint into 
the partition function [I] . This entangled polymer system 
has a similar mathematical structure to the trajectory of 
a classical particle under a magnetic field following the 
Biot-Savart law, which allows us to analyze the polymer 
entanglement in favor of the magnetic-field induction un¬ 
der the constraint of integer linking number. This leads 
us to construct BE theory [31] that avoids self-linking 
contrary to Chern-Simons theory. 

We agree that it is extremely difficult to solve such 
complicate equations obtained from Chern-Simons the¬ 
ory exactly and provide the correspondence of parame¬ 
ters between the microscopic and phenomenological the¬ 
ories 0 [11]. In this respect we adopt a new strat¬ 
egy to make a link between the phenomenological tube 
model and the microscopic BE theory. Inspired by previ¬ 
ous works HHH], we derive an effective Maxwell-London 




equation for the dynamics of entangled polymer melts, 
starting from the topological BF theory. We interpret 
the characteristic features of the tube model with this 
Maxwell-London equation, which leads us to identify the 
dynamics of entangled polymers with the Higgs phenom¬ 
ena. 


Introducing an artificial electromagnetic vector potential 
A/(r) given by other polymers except for the 


V X A 7 (r) = ^jj(r), (3) 



where the most general solution to the Eq. © is 



(4) 


it is convenient to express the winding number of the 
polymer with others as follows 


Nj = f dr'^Aj{r) ■ jj{r). (5) 

I 

If we consider only two polymer chains for simplicity, it 
is straightforward to see how Eq. ([^ shows the winding 
number 


FIG. 1: Schematic illustration of topological constraint of 
polymer melts. We assume that the number of entanglement 
is preserved during time scale we are interested in. 


iV2 = y drA2(r) • j2(r) 


Topological field theory: The polymer density and cur¬ 
rent for N chains of length L are given by 


N i-L 

P(r) = E / 

jw=i:/ 


ds6{xi{s) — r) 


dsx{s)6{xi{s) — 


( 1 ) 

r). (2) 



X ■ h{r)dr 


dr'{s) X dr”(s') ■ 


r'js) - r”{s') 
|r'(s) — r"(s')|3 


( 6 ) 


Introducing this topological constraint into the partition 
function for the polymer, we obtain 


Z7 = 


oo oo N 

-fiiNi 


Ni=0 


E <^(E - Ni) ( d V X A 7 (r) - ^ j,(r) A 7 - / dr^, Aj{r) ■ jj{r) 


N,i=0 j^I 
OO N 


= E^ 


-fiiNi 

Ni=0 Nij=0 j^I 


Nil 




j UfL,DX,{s)DAi{r)DCi{r) 


exp{-(/ ds'^^X]{s)-igi J d‘^rj/(r) • A/(r) - f J d^rCi{r) ■'^ij{r) + i j d‘'rC/(r) • [V x A 7 (r)]) |, 
2=1 

(7) 


N 




where gi is a Lagrange multiplier field to impose the 


constraint SiNi — A 7 (r) • j 7 (r) 


and C 7 (r) 



























System 

Superconductor (SC) 

Entangled polymer 

Medium 

SC vacuum 

Entangled polymer complex 

SC boundary 

Local averaged position of a chain 

Degrees of freedom 

Current 

Tangential flux of a chain 

External magnetic field 

Test chain 

Magnetic field induction 

Flux of entangled chains 

Topology constraint 

Flux quantization 

Preserving linking number 

Physics 

Photon mass 

Retracting force at intersections 

Magnetic field penetration 

Diffusion of tangential flux 

Decay length 

Tube radius 


TABLE I: An analogy between polymer entanglement (tube model) and superconductor 


is an auxiliary gauge field to impose the constraint 
i5 X A/(r) — integer to repre¬ 

sent a Gaussian linking number between the polymer 
and jth one (j ^ I). Nj = Nij is an integer to 

express a total Gaussian linking number between the 
polymer and all others (j 7 ^ /). fij is an energy cost per 


a linking event. Although Gaussian-fluctuation polymers 
have been assumed, the formal expression can be gener¬ 
alized to take their interactions. 


It is straightforward to write down the canonical parti¬ 
tion function for N polymers, given by 


Zm = 


00 00 N tjN at \ ^ ^ TT^ Af I /‘oo 

Ng=0 Ni= 0 1=1 Nij=0 j^I 


Nij—0 

N ,L 


+z J d'^rC/(r) ■ [V X A/(r)]^|, 


7=1 


N 


U^^^DXj{s)DAi{r)DCj{v)expl-^( I ds^Xj{s) - ^ ( d^r5/j/(r) • A/(r) - z f (f^rC/(r) • ^ (r) 


0 




( 8 ) 


where Nq = i® ^^e total Gaussian linking num¬ 

ber and hg is an energy cost per a linking event. The 
partition function may be more manageable in a form of 
the grand canonical partition function, given by 


Zgc = 

Ar=o 




(9) 


where /r is a chemical potential for a polymer. 

An essential feature in this formulation lies in the so 
called BF term, given by 


which describes mutual entanglement of polymers. In 
addition to the BF term, the dynamics of the fluctuating 
polymers will give rise to the Maxwell dynamics for the 
emergent gauge fields of A/(r) and C/(r). This may be 
regarded as the effect of renormalization for the dynamics 
of photons from the dynamics of electrons in the theory of 
quantum electrodynamics m- As a result, we obtain the 
so called BF-Maxwell theory as an effective description 
for the dynamics of entangled polymers, given by 


N 


Sbf = i d'^vY,Ci{v)-[VxAi{v)l (10) 


7=1 




J -B^^(r)] -b (r)] -b *C/(r) • [V x A/(r)]^, 


( 11 ) 






















where both effective coupling constants of and for C/(r) and 

are proportional to gj since these contributions originate 
from the fluctuating polymers. Considering 


^ , d AdX,(s) d 

^^ 7 -= I 


^ ds dXJs) 
1—1 ' ' 

J ds^Jd‘^rS^^\r-X,is))±,is)-Vr 


ds 

L N 


i=l ■ 


we obtain 


EA^(r) = -dtAj(r) = 
BAj(r) = V X A/(r), 


N 

^j,(r) • VA/(r), 
2=1 


( 13 ) 


N 

z=l 

Then, effective electric and magnetic fields are given by 


for A/(r). 


Ec,(r) = -atC,(r) = 
Bc,(r) = V X C,(r) 


N 

2 = 1 


( 12 ) 


Maxwell-London equation: Applying the least action 
principle to the BF-Maxwell action, it is straightforward 
to derive the following equation of motion 


N N 


4 E E[j*W ■ W • ^]C/(r) + ^y^Ciir) - 4 vv • C,(r) = 


2=1 2' = 1 
N N 


^ ^[j,(r) • V][j,,(r) . V]A,(r) + 4 v'A,(r) - ^VV • A,(r) = 


-iV X A/(r), 
-iV X C/(r). 


i=i i'—i 


( 14 ) 


If we focus on the linear regime for gauge dynamics, we 
are allowed to neglect the electric-field term. The reason 
is that the polymer current can be expressed in terms 
of gauge potentials, referred to as constituent equations 
and giving rise to higher-order dynamics of gauge fields. 
This may be regarded to be nothing but the Ohm’s law 
in metals although it is not clear how to construct precise 
constituent equations in the present situation, where the 
dynamic information of the entangled polymers should 
be incorporated. 

Taking the Coulomb gauge of V • C/(r) = 0 and V • 
A/(r) = 0 , we find the Maxwell-London equation 

4v2Bc,(r) = e2Bc,(r), 

4v2BA,(r) = e2BA,(r). (15) 

Ca 

Recall 

N 

BA^(r) = EJj W’ Bc,(r) = (16) 

We write down the Maxwell-London equation in a more 


snggestive fashion 

(DV^ - M 2 (r)) B(r) = 0 , ( 17 ) 

where the diffusion coefficient of D is introduced for phys¬ 
ical interpretation. It is interesting to notice that M(r) 
can be identified with an effective mass of photons in¬ 
side SC that corresponds to an inverse of the penetra¬ 
tion depth of magnetic fields. In the setup of the BF 
theory it is proportional to gf, where gj is the coupling 
function to enforce the constraint of the Gaussian link¬ 
ing number. Within our semiclassical framework, it can 
be determined by the saddle-point approximation for the 
resulting free energy. This procedure may be one of the 
most complex parts in our analysis. In the present study 
we do not perform this saddle point analysis and take the 
mass parameter as a function of the distance away from 
the average position of a test polymer. In the language of 
polymers, the penetration depth is read as a decay length 
of the transverse distribution of the test polymer current 
from its local averaged position. The mass may be under¬ 
stood as a retracting force at the intersections along the 
contour. In SC, the dissipationless current is induced to 
cancel the magnetic field inside as seen in Fig.j^ leaving 
the shallow penetration depth at the surface. In polymer 
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FIG. 2: The transverse distribution of the chain is plotted at 
very low polymer density of 0.01. The distribution is fitted 
with Gaussian perfectly except the center of local average 
position where the attraction of Lennard-Jones interaction 
affects. 




Np Np 

FIG. 3: In (a), (b) the transverse distribution of the chain 
from the primitive path at different Np, and fcg. In (c), (d) 
the areal fractions of Gaussian distribtuion and the exponen¬ 
tial decaying function are plotted as function of Np when the 
transverse distributions are htted with linear combination of 
Gaussian distribution and exponential decaying function. 


melts, entangled chains (with a test chain) restrict large 
deviations of the test chain segment in the normal direc¬ 
tion. Thus, we conclude that the decay of the transverse 
polymer current in the normal direction is an essential 
character resulting from the topological constraint. This 
picture is well incorporated in the tube model. In com¬ 
parison with the tube model, the backbone of the tube 
is a local average position, and the radius of the tube 
corresponds to the decay length (the penetration depth 
in SC). In this argument, the exponential decay of the 
chain leakage out of the tube is well understood. 

Comparison with numerical results 

Method: The idea, that the transverse chain flux dis¬ 
tribution is a linear combination of exponential decay¬ 
ing function and the Gaussian distribution, is tested by 
computer simulations. The linear potential to the nor¬ 
mal direction of local average configuration will cause the 
exponenital decay of transverse current flux. A compu¬ 
tational time of the regular MD or MC simulations for 
entangled long polymer melts takes very long. Consider¬ 
ing the relaxation time, which scales as about ( Np 
is polymerization degree.), we only are able to simulate 
short chain of polymer {Np < 1000). 

Under this restriction, we performed molecular dynamics 
simulations based on the coarse-graining entanglement 
model developed by Kremer and Crest [55] . Temperature 
of the system is kept to be 297K using the Langevin 
equation. 

Three interactions are considered to reflect the monomer 
size, chain connectivity, and the bending rigidity of the 
polymer. The interaction between monomers is incorpo¬ 


rated in a Lennard-Jones (LJ) potential. 





r < Tc 
r >rc 


(18) 


where, Tc = 2 5 ctq . ctq is a length unit in LJ model repre¬ 
senting a diameter of monomer. 

The chain is sustained by a FENE potential. 


Ufene 


-QNkRl In (1 - (r/i?o)2) 
oo 


r < Rq, 
r > i?o 


(19) 


where i?o = 1.5 (To. e is an energy unit in LJ potential. 
To adjust the stiffness of the chain, we add a bending 
potential. 


Ubend = ke{l- cos 9) (20) 

where, ke is the bending stiffness, and the angle 9 is de¬ 
fined as cos 9i = . When kg = 0, the 

chain is flexible. Chain becomes stiff increasing the bend¬ 
ing stiffness, for example at kg = 2, the characerisitic 
ratio Coo is about 3.4 considered as semiflexible poly¬ 
mer. According to the reference [55] N^ is 65 for kg = 0, 
and 23 for kg = 2. As Kremer and Crest pointed out 
at e = ksT, and k = 30-%, there is a huge energy bar- 

rier (~ lOksT) for two chains to cross. This ensures the 
uncrossbility of the chain. 

The time step is chosen as tq = 0.01. The friction coeffi¬ 
cient for Gaussian noise is set O.bm/ro. To properly equi¬ 
librate the system, we adjust the simulation time scale 
to exceed the longest relaxation time (the disengagement 
time, Td = 4.5(A^^/A^e)T'o)- 10® time steps are used for 






















equilibration, and another 10® time steps are used for the 
data production. 

The transverse distribution of the chain from its local av¬ 
eraged position is measured within the constraint release 
time. For this purpose, we apply the direct tube sam¬ 
pling (DTS) [20]. At every 5 x 10® steps (which exceeds 
Rouse time, Tp> ~ l.SiV^ro) during production stage, we 
write out the conformation of polymer. In DTS, we find 
the primitive paths of each conformation are saved during 
production stage. In order to keep the topological con¬ 
straint preventing the reptation motion, disengagement 
or constraint release, we fix the ends of the polymers. 
Then, by cooling the temperature of the system to zero, 
we find out the energy minimized conformation of the 
polymer which is called primitive path. Then, we return 
back to the normal temperature and measure the devia¬ 
tion of the chain from the primitive path. 

Results: Linking number is preserved only for an inter¬ 
mediate time scale ^ when the topological 

constraint affects the chain dynamics but the reptational 
diffusion, disengagement or constraint release are still 
suppressed. Incorporating the dynamical aspects of the 
entangled polymers into the static Maxwell-London equa¬ 
tion for this regime, we propose an extended Maxwell- 
London equation 

dtB{r,t) = - M^ir)B{r,t), (21) 

where 

B(r,f) oc j(r,t) = J ds-^X{s,t)6^^'>(X{s,t) - r). 

( 22 ) 

Within this time scale, the distribution of the intersec¬ 
tion is not fully uniform over the chain, but localized. For 
the part of the chain far from the intersection, the diffu¬ 
sive motion would be dominant, i.e. the polymer current 
of Ji will follow the Gaussian distribution, described W 
dtB{r,t) = I?V^B(r,t) especially, near |r| « 0. In Fig|R 
this idea is verihed from the numerical simulation. At 
very low polymer density (p = 0.01), the transverse dis¬ 
tribution fits perfectly with a Gaussian function except 
at the vicinity of the local average position where the 
attraction of Lennard-Jones localizes the monomer dis¬ 
tribution. The diffusion time scale (or relaxation time 
at quenched intersections) is proportional to the chain 
length per the entanglement {td — V^<^o) which is 
nothing but width of the Gaussian distribution. On the 
other hand, near the intersections the fluctuations of the 
test polymer will be suppresed. As a result, we expect a 
static solution in which the mass term is dominant, de¬ 
scribed by the previous Maxwell-London equation. The 
decay length Ad is about In a nutshell, the tube 

diameter can be justified as a dominant decaying length 
scale between rjj and A^. 


In our derivation, the mass term is proportional to the 
energy of linking events, and its spatial frequency (linking 
number), i.e. the mass term is dominant when the inter¬ 
action energy at intersection is strong and linking num¬ 
ber is large. One simple way to achieve both conditions 
is turning the bending energy up. In Fig. |^(a),(b) we 
plot the transverse distribution of the chain as a normal 
distance a from its local averaged conhguration varying 
ke (at p = O.SScTf)"®). As we noted, increasing fcg, a chain 
transits from a flexible chain to a semi-flexible chain. In 
accordance, the linking number and interaction energy 
at the intersection increases. Indeed the transverse dis¬ 
tribution of current decays faster for kg = 2 regardless 
Np (Fig.§(a),(b)). 

For each distribution, we fit them with a linear combi¬ 
nation of the Gaussian distribution and the exponential 
decaying function. They are fitted very well with this 
function. Then, the areas of the Gaussian part and the 
exponential decaying part are plotted in Fig. §(c),(d). 
This shows which contribution is dominant to the distri¬ 
bution. At shorter polymer length {Np = 75,200), the 
exponential decaying part is dominant regardless of kg. 
This is because of the artificial retracting force from the 
quenched end points. The artificial effect from end points 
gets weaker at larger chain length. In polymer lengths 
of Np = 300 ~ 1000, the trends are clearly splitted as 
kg. At kg = 0 (flexible polymer with larger A’g), the 
Gaussian part is dominant. The trend is a bit dependent 
on the chain length, and the Gaussian part is slightly 
more important for longer chains. This can be under¬ 
stood from the weaker contribution from the fixed ends 
for longer chains. In contrast, the exponential decaying 
part is more important for kg = 2. The contribution is 
larger for longer chains. This is because longer chains 
can form tighter networks, and thus the interaction en¬ 
ergy at intersections is stronger. Indeed as we expected, 
the interaction energy at the intersections and the linking 
number increase the mass term. 

Towards the second “quantization”: Our phenomenolog¬ 
ical Maxwell-London equation based on the topological 
BF theory seems to be not inconsistent with the suc¬ 
cessful tube model. However, there still exist important 
unanswered questions, which may be difficult to be solved 
within the present formulation. For example, it is not 
clear how to express the mass parameter as a function of 
measurable physical quantities. We point out that this 
fundamental difficulty of our formulation results from the 
fact that the theory is written in the language of the 
first “quantization”. We need to represent the BF theory 
in the second-quantization formulation, where not only 
gauge dynamics but also polymer dynamics is expressed 
in terms of field variables instead of the position vari¬ 
able. This work has been performed before |29|, where 
the generating function of 





leads to the Edwards-Anderson equation 


Here, the s) field satisfies 


dsq{r, s) = BV^q(r, s) (24) 

in the saddle-point analysis and an internal energy 


s) = -BV^q^r,s). 


(29) 


E = - lim — In W. (25) 

N—^OO N 

The density and current of polymers are given by 
p(r,s) = g'l'(r,s)g(r,s) 

and 

j(i’,s) = -^(' 7 ^(i-,s)[Vrg(r,s)] - [VrQ^(r, s)]g(r, s)), 

(27) 


Since the equation of motion for the q{r, s) field differs 
from that of the q^ {r, s) field, it is convenient to introduce 
^ 25 ) a Nambu-spinor [30] 


■!/>(r,s) = 


g(r,s) 

qHr,s) 


(30) 


respectively, which satisfy the conservation law 

Then, the generating function can be reformulated as 
dsp{r, s) + Vr • j(r, s) = 0. (28) follows 


W = 


£)V’^(r, s)£>V'(r, s)exp|-y ds J s){dsT3 - D\7l)tp{r, s)Y 


(31) 


where the density and current of polymers are expressed respectively, 
by 

P(i’,s) = iV'’^(r,s)V’(r,s) (32) 


and 

j(i’,s) 


1 

Id 


(^'(/''''(r,s)T3[VrV'(r,s)] - [Vr'0'''(r,s)]T3'0(r, 


Introducing the chemical potential (/r) and effective in- 
/ teractions (V) of polymers, we can describe interacting 
(ipi)lymers at finite density, given by 


W = 


D7/!'^(r,s)D'ip(r,s)exp - J ds J s){dsT 3 - p - D\7l)ip{r, s) + s)tlj{r, .(34) 


Entanglement of polymers can be formulated as before, introducing gauge fields as follows 









(35) 



-K^^ip^{r,s)tlj{r,s)[Ar{r,s)]'^ + 


Here, k is the Lagrange multiplier field to impose the 
Gaussian linking number and Nq is the total number 
of the Gaussian linking number. In this formulation we 
allow self-linking of a polymer in addition to mutual en¬ 
tanglement, described by the Ghern-Simons term 

Scs = J ds J 

(36) 

where 0 is a statistical angle, which can be used as a 
phenomenological parameter. By substituting the time 
derivative dt into the segment derivative ds, we can 
revisit the previous effective Maxwell-London equation 
from the effective Ghern-Simons field theory. 

This field-theoretic formulation allows us to introduce 
renormalization effects of both dynamics of gauge and 
polymer fields self-consistently at least in the level of 
random phase approximation (RPA) [3T] . Although this 
Ghern-Simons field theory has a similar structure with 
that of fractional quantum Hall effect, we can’t apply 
its solution directly. There exists a substantial differ¬ 
ence between two Ghern-Simons field theories: The time 
derivative does not have an i factor, where f is a complex 
number with = —1 |32j . Self-consistent renormaliza¬ 
tions within the RPA level are required near future. 
Discussion: Starting from the topological BF theory in 
the first-quantization representation, we derive the effec¬ 
tive Maxwell-London equation. We connect the solution 
of this phenomenological equation with the physics of 
the tube model. The essence of the tube model is ex¬ 
plained by the Meissner effect of our phenomenological 
equation. We revisit this formula with an alternative ap¬ 
proach based on second quantization. These results seem 
to be consistent with our numerical analysis. 

We would like to emphasize that the role of the topolog¬ 
ical constraint (the preserved Gaussian linking number) 
differs from that of effective interactions between seg¬ 
ments of polymers. Effective interactions between seg¬ 
ments of polymers can be translated into dynamics of 


D^^{r, s)Dip{r, s) 

(r, s)r 3 ]T 3 - fj.- D[Vr - f«:A(r, s)r3]^) V’(i', «) 

0 1 ' 

)V’(r, s)]2 -f i — e^,,xAf,{r, s)a^AA(r, s)| . 


longitudinal gauge fluctuations, i.e., Coulomb interac¬ 
tions, which would be screened by polymer fluctuations, 
nothing but the Debye screening in metals |33j and thus, 
allowed to be neglected at low energies. This massive 
dynamics of longitudinal gauge fluctuations should be 
distinguished from Higgs phenomena, where transverse 
gauge fluctuations become massive, realized in the SC 
media. The phenomenological Maxwell-London equation 
implies that the tube model should be understood within 
the presence of the topological constraint beyond effec¬ 
tive interactions between polymer segments. 

The Ghern-Simons field theory allows us to pursue an 
analogy between superconductivity and entangled poly¬ 
mer complex in depth. The fundamental concept of U(l) 
symmetry breaking in superconductivity implies the gap¬ 
less Goldstone mode and gapped Higgs mode, where the 
gapless sound mode is pushed up to the plasmon mode 
when there exist transverse gauge fluctuations [33]. In 
order to search such plasmon and Higgs modes in the 
entangled polymer complex, we need to investigate the 
dynamics of entangled polymers, integrating over gauge 
fluctuations. This deep connection can be realized in the 
second-quantization representation. Then, the dynamics 
of corresponding the plasmon and Higgs modes would 
reveal characteristic responses of the entangled polymer 
complex such as compressibility, viscosity, and etc. These 
will be our future direction. 
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